Universal Selftrapping in Nonlinear Tight-binding Lattices 
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We show that nonlinear tight-binding lattices of different geometries and dimensionalities, dis- 
play an universal selftrapping behavior. First, we consider the single nonlinear impurity problem 
in various tight-binding lattices, and use the Green's function formalism for an exact calcula- 
tion of the minimum nonlinearity strength to form a stationary bound state. For all lattices, 
we find that this critical nonlinearity parameter (scaled by the energy of the bound state), in 
terms of the nonlinearity exponent, falls inside a narrow band, which converges to e ' 2 at large 
exponent values. Then, we use the Discrete Nonlinear Schrodinger (DNLS) equation to examine 
the selftrapping dynamics of a single excitation, initially localized on the single nonlinear site, 
and compute the critical nonlinearity parameter for abrupt dynamical selftrapping. For a given 
nonlinearity exponent, this critical nonlinearity, properly scaled, is found to be nearly the same 
for all lattices. Same results are obtained when generalizing to completely nonlinear lattices, 
suggesting an underlying selftrapping universality behavior for all nonlinear (even disordered) 
tight-binding lattices described by DNLS. 



The Discrete Nonlinear Schrodinger (DNLS) equation 
is a paradigmatic equation describing among others, dy- 
namics of polarons in deformable mediaQ, local modes 
in molecular systems and power exchange among non- 
linear coherent couplers in nonlinear optics [B. Its most 
striking feature is the possibility of "selftrapping" , that 
is, the clustering of vibrational energy or electronic prob- 
ability or electromagnetic energy in a small region of 
space. In a condensed matter context, the DNLS equa- 
tion has the form 
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where G n is the probability amplitude of finding the elec- 
tron (or excitation) on site n of a d-dimensional lattice, 
e n is the on-site energy, V is the transfer matrix element, 
Xn is the nonlinearity parameter at site n and a is the 
nonlinearity exponent. The prime in the sum in (m) re- 
stricts the summation to nearest-neighbors only. In the 
conventional DNLS case, a = 2 and Xn is proportional to 
the square of the electron-phonon coupling at site n. Q 

Considerable work has been carried out in recent years 
to understand the stationary and dynamical properties 
of Eq. (II) in various cases. In particular, we point out 
the studies on the stability of the stationary solutions 
in one and two dimensions for the homogeneous case 
(e n = 0,Xn = x)f| [§, the effect of point linear impu- 
rities on the stability of the 2-D DNLS solitonsM, the 
effects of nonlinear disorder (e n = 0, Xn random)|p| and 
of linear disorder (x„ = x, e n random) || on the selftrap- 
ping dynamics of initially localized and extended exci- 
tations in a chain. The results obtained in these studies 
suggest that, in general, the effect of nonlinearity is quite 
local for initially localized excitations, and that disorder 
leaves the narrow selftrapped excitations unaffected, al- 
though it does affect the propagation of the untrapped 



portion ("radiation"). In this Letter we show that, for 
an initially localized excitation, the dynamics of selftrap- 
ping in various different lattices of different dimensional- 
ities, is universal and depends mainly on the nonlinear- 
ity strength at the initial site, the nonlinearity exponent 
and the coordination number, and much more weakly on 
other topological features of the lattice. 

Bound states. A tight correlation has been observed be- 
tween the existence of bound states for a given nonlinear 
lattice and the ability of the lattice to selftrap an initially 
completely-localized excitation: the critical nonlinearity 
strength for dynamical selftrapping is always greater than 
the one needed to produce bound state(s). We begin by 
showing that the minimum nonlinearity needed to pro- 
duce a bound state in different lattices, shows universal 
features. 

We consider the problem of determining the bound 
state for an electron in a d-dimensional homogeneous lat- 
tice that contains a single generalized nonlinear impurity 
at the origin n = 0. The Hamiltonian is H — Hq + &i> 
where H = V J2 n «( l n )( m l + n - c - ) * s * ne unperturbed 
tight-binding Hamiltonian with hopping constant V and 
Hi = x\Co\ a 1 0) (0 1 corresponds to the nonlinear im- 
purity perturbation. The {|n}} represent Wannier elec- 
tronic states, and we have set e n = 0. For convenience 
we normalize all energies to a half bandwidth, B and 
define: z = E/B, H = H/B and 7 = X /B. The di- 
mensionless lattice Green function G = l/(z — H) can 
be formally expanded as|§ G = G (0 ) + G^HiG^ + 
G (0) i/iG (0) iJiG(°) + ..., where G<°) is the unperturbed 
(7 = 0) Green function and H x = 7|G | a |0>< 0|. The 
sum can be carried out exactly to yield 
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where G mn = (m|G|n). The energy of the bound 
state(s), Zb is obtained from the poles of G mn , i.e., 



by solving 1 
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The bound state am- 



plitudes G„ are obtained from the residues of G mn (z) 



at z = Zb- In particular, |C 
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^oo ( z b) / G qq' (zb). Inserting this into the bound state 



energy equation leads to 
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We proceed to solve (tfl) numerically, using the exact, 
known expressions for G 00 for several lattices|lQ, pT|: 
one-dimensional (1-D), square, triangular, simple cubic 
and Bethe lattices with connectivities 3, 5 and 100. This 
allows us to compare lattices with different dimensional- 
ity, coordination number Z, length of shortest loops, etc. 
In general, for a given a value there will be a minimum 
value of x below (above) which, there is (are) no (two) 
bound state(s). Just at the critical nonlinearity value, 
we obtain exactly one bound state. The exception is the 
1-D lattice where one needs in addition, a > 2 1131 . 




FIG. 1: The critical nonlinear paremeter, "f c /zb, for 
Bound States. Thick lines correspond to Cubic, Square 
and ID cases. Thin lines correspond to Bethe lattices 
with K — 3, 5 and 100 in ascending order near a = 0. 
Dotted lines represent the Triangular cases: sgn(x/V) > 
(upper line) and sgn{x/V) < 



Figure 1 shows the critical nonlinearity parameter 7 C , 
scaled by the energy of the bound state, in terms of 
a, the nonlinearity exponent, for all the lattices exam- 
ined. These curves are independent of sgn(x/V), ex- 
cept for the triangular lattice, due to the asymmetry 
of its Green function with respect to the energy vari- 
able. In this case there are two curves depending on 



sgn(x/V). All curves in Fig.l fall inside a "band" which 
narrows as a increases, converging towards a constant 
value. To calculate it, we solve (0) exactly in two cases: 
the one-dimensional lattice [g2| and the Bethe lattice in 
the limit of infinite connectivity (numerically indistin- 
guishable from K = 100). In both cases we obtain: 



lim ( — 

Zb 



= 1/2 



1.65. 



(4) 



We have traced the validity of (|j) for the other lattices 
up to high a values (10 3 for the square and cubic lattices; 
10 5 for the rest) with no discernible deviation. 



Selftrapping Dynamics. We now examine the ability of a 
given lattice to dynamically selftrap an excitation, orig- 
inally placed completely on the impurity site, by com- 
puting the minimum nonlinearity needed to give rise to 
abrupt selfptrapping. The time evolution is given by 
Eq.([j]) with e n — and Xn — X &i,o- The numeri- 
cal scheme is that of a fourth-order Runge-Kutta, where 
the accuracy is monitored through total probability con- 
servation. To avoid undesired boundary effects, a self- 
expanding lattice is usedB. To ascertain the presence or 
absence of a sharp selftrapping transition, we compute 
the long-time average probability at the impurity site, 
defined by 



P = lim (1/T) / 

T^oo V J 



\C Q (t)\ 2 dt, 



|q,(0)| = 1. (5) 



Typically, Pq vanishes for nonlinearity parameters be- 
low a critical value Xc and the particle escapes from the 
impurity site in a ballistic manner. This is determined 
from an examination of the excitation's mean square dis- 
placement (u(t)) = Yin n2 IOi| 2 - For nonlinearity val- 
ues greater than Xc Pq remains finite and increases with 
X, converging towards unity at large x- The untrapped 
portion escapes to infinity, also in a ballistic manner, but 
with a much lower "speed" y/(u(i))/Vt. Thus, from the 
examination of Pq we determine the critical nonlinear- 
ity parameter Xc for dynamical selftrapping (usually for 
Po«l/2). 

For a particular lattice and a given exponent a, we nu- 
merically determine the critical nonlinear parameter Xc-, 
scaled by Eb (where Eb is the unnormalized bound state 
energy correspondig to this Xc) for abrupt selftrapping. 
Figure 2 shows Xc/Eb for all the lattices examined, and 
for several a values that give rise to sharp selftrapping 
(for a < 1, the selftrapping is not sharp). We see that, 
for the wide range of geometries and dimensionalities in- 
volved, this critical (dynamical) nonlinearity is nearly 
independent of the lattice and increases monotonically 
with the nonlinearity exponent. This is specially true in 
the all-important conventional DNLS case (a — 2). It 
would seem that, in the a regime where abrupt selftrap- 
ping takes place (a > 1), the only relevant parameters 



are the nonlinearity at the impurity site and the coordi- 
nation number of the lattice. The rest of the topological 
features is of secondary importance. In all cases, with the 
exception of the triangular lattice, the critical nonlinear- 
ity is independent of the sign of x/V ' . For the triangular 
lattice we note that Xc/Et gets shifted a bit upon chang- 
ing the sign of x/V ■ This probably trails back to the 
asymmetry of the unperturbed triangular lattice's Green 
function Gq under a sign change of its argument pj . 
All the rest of the lattices are symmetric in that respect. 
The increase of \ c with a is to be expected since, in 
the continuum limit , increasing a is equivalent to in- 
creasing the dimensionality of the system || [b|; this in 
turns increases the effective coordination number mak- 
ing it harder to selftrap the excitation; hence, the need 
for larger nonlinearities. Also, the obtained values of Xc 
in the dynamical case are all higher than for the bound 
state case, confirming the conjecture that the onset of the 
stationary bound state is a precursor for dynamical self- 
trapping. However, the lack of a superposition principle, 
makes it hard to establish formally the (observed) con- 
nection between the dynamical and the stationary DNLS 
problem. An alternative normalization for Xc is to use 
the half bandwidth B instead of E\,. In that case, all the 
curves in figure 2 lose a bit of flatness, but the tendency 
is otherwise unaltered. 
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Figure 3 shows the results obtained for the critical nonlin- 
earities. The curves are virtually the same as the ones in 
Fig. 2. (The case a = 1 does not display abrupt selftrap- 
ping like the rest, thus Xc is not precisely defined here). 
This is due to the fact that, once the abrupt selftrap- 
ping is set, most of the probability is on the initial site, 
which gives, by conservation of probability, very small 
probability amplitudes for the rest of the lattice sites, 
making their nonlinear contribution negligible: they have 
become, in fact, linear for all selftrapping purposes and, 
in this way we are back to the single nonlinear impurity 
results. The greater the a value, the closer the system 
to the nonlinear impurity case. This is vividly illustrated 
by the limiting curves for large a in Figs. 2 and 3., which 
coincide. 
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FIG. 3: The dynamic critical nonlinearity parameter 
scaled by the bound state energy (from the stationary 
impurity problem) for completely nonlinear lattices. The 
values for the nonlinearity exponent are a = 1,2,3,4 
from bottom to top (the case a = 1 does not show abrupt 
selftrapping; hollow symbols represent sgn(x/V) < 0). 
The limiting (upper) curve correspond to a — 1000. 



FIG. 2: The dynamic critical nonlinearity parameter Xc 
scaled by the bound state energy for one nonlinear im- 
purity in various lattices. The values for the nonlinearity 
exponent are a — 1,2,3,4 from bottom to top (hollow 
symbols represent the case sgn(x/V) < 0). The limiting 
(upper) curve correspond to a = 1000. 



This characterization of the selftrapping properties of 
nonlinear tight-binding lattices of different geometries 
and dimensionalities, in terms of a single parameter, 
namely the bound state energy for the one-impurity 
problem (or the half bandwidth B for quick estimations) , 
could be useful in several areas, given the paradigmatic 
character of DNLS. 
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